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Abstract. In this paper we develop new applications of variational analysis and generalized differ- 
^v^j ■ entiation to the following optimization problem and its specifications: given n closed subsets of a 

Banach space, find such a point for which the sum of its distances to these sets is minimal. This prob- 
lem can be viewed as an extension of the celebrated Fermat-Torricelli problem: given three points on 
the plane, find another point such that the sum of its distances to the designated points is minimal. 
The generalized Fermat-Torricelli problem formulated and studied in this paper is of undoubted 
mathematical interest and is promising for various applications including those frequently arising in 
location science, optimal networks, etc. Based on advanced tools and recent results of variational 
analysis and generalized differentiation, we derive necessary as well as necessary and sufficient op- 
timality conditions for the extended version of the Fermat-Torricelli problem under consideration, 
which allow us to completely solve it in some important settings. Furthermore, we develop and 
justify a numerical algorithm of the subgradient type to find optimal solutions in convex settings 
and provide its numerical implementations. 
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1 Introduction and Problem Formulation 

In the early 17th century Pierre de Fermat proposed the following problem: given three 
points on the plane, find a fourth point such that the sum of its Euclidean distances to the 
three given points is minimal. This problem was solved by Evangelista Torricelli and was 



named the Fermat-Torricelli problem. Torricelli's solution states the following: if none of 
the interior angles of the triangle formed by the three fixed points reaches or exceeds 120°, 
the minimizing point in question is located inside this triangle in such a way that each side 
of the triangle is seen at an angle of 120°; otherwise it is the obtuse vertex of the triangle. 
This point is often called the Fermat-Torricelli point. 

In the 19th century Jakob Steiner examined this problem in further depth and extended 
it to include a finitely many points on the plane. A number of other extensions have been 
proposed and studied over the years. This and related topics have nowadays attracted 
strong attention of many mathematicians and applied scientists; see, e.g., [3, 4, 11, 20] with 
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the references therein for the history, various extensions, modifications, and applications 
to location science, statistics, optimal networks, etc. Note that, despite beautiful solutions 
obtained for particular extensions of the Fermat-Torricelli problem, we are not familiar 
with theoretical methods and/or numerical algorithms developed in rather general settings. 
We particularly refer the reader to [6, 10, 21] and the bibliographies therein to Weiszfeld's 
algorithm and its modifications for the problem of minimizing weighted sums of Euclidean 
norms (also known as the Weber problem) and to [1] for efficient interior point type methods 
for similar problems in finite-dimensional spaces. 

In this paper we study a far-going generalization of the Fermat-Torricelli problem that 
is formulated below. It extends, in particular, a generalized version of the classical Steiner 
(and Weber) versions with replacing given points therein by a finitely many closed sets in 
Banach spaces. Furthermore, our new extension of the Fermat-Torricelli problem covers 
a fast majority of the previous ones and seems to be interesting for both the theory and 
applications to various location models, optimal networks, wireless communications, etc. 

We propose to employ powerful tools of modern variational analysis and generalized 
differentiation to study the extended version of Fermat-Torricelli problem and its specifi- 
cation from the theoretical/qualitative and numerical/algorithmic viewpoints. In the first 
direction our goal is to derive necessary as well as necessary and sufficient optimality condi- 
tions for generalized Fermat-Torricelli points and then to use them for explicit determining 
these points in some remarkable settings. Our numerical analysis involves developing an 
algorithm of the subgradient type and considering its specifications and implementations in 
the case of the generalized Fermat-Torricelli problem determined by an arbitrary number 
of convex sets in finite-dimensional spaces. 

Let us now formulate the generalized Fermat-Torricelli problem of our study. Consider 
the so-called minimal time function 

Tg(x) := inf {t > 0| On (x + tF) / 0} (1.1) 

with the constant dynamics x G F described by a closed, bounded, and convex subset F ^ 
of a Banach space X and with the closed target set U ^ in X; these are our standing 
assumptions in this paper. We refer the reader to [7, 8, 9, 15] and the bibliographies therein 
for various results on minimal time functions and their applications. When F is the closed 
unit ball IB of X, the minimal time function (1.1) becomes the standard distance function 

d(x;fi) =inf{||x-w|| | co G Q} (1.2) 

generated by the norm || • || on X. Given now an arbitrary number of closed subsets / 0, 
i = 1, . . . , n, of X, we introduce the generalized Fermat-Torricelli problem as follows: 

n 

minimize T(x) := J^T(f.(x), x € X. (1.3) 
i=l 

For F = IB in (1.3) this problem reduces to 

n 

minimize D(x) := d(x; fij), (1-4) 
i=i 
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which corresponds to the Steiner-type extension of the Fermat-Torricelli problem in Banach 
spaces when all the sets fij, i = 1, . . . , n, are singletons. Observe that even in the latter 
classical case the optimization problem (1.3) and its specification (1.4) are nonsmooth while 
being convex if all the sets Jlj have this property. It is thus natural to study these problems 
by means of advanced tools of variational analysis and generalized differentiation. 

The rest of the paper is organized as follows. In Section 2 we define and discuss basic 
tools of variational analysis needed for formulations and proofs of the main results of this 
paper. Section 3 is devoted to computing and estimating subdifferentials of minimal time 
functions that play a crucial role in our study of the generalized Fermat-Torricelli problem 
and its specifications. In Section 4 we derive necessary conditions for the general problem 
(1.3) in Banach spaces as well as necessary and sufficient conditions in the case of convexity, 
which are then used for complete descriptions of Fermat-Torricelli points in some important 
settings. Finally, Section 5 presents and justifies a numerical algorithm for solving the 
generalized Fermat-Torricelli problem with convex data in finite-dimensional spaces. 

Throughout the paper we use standard notation and terminology of variational analysis; 
see, e.g., [13, 18]. Recall that, given a set-valued mapping G: X =f X* between a Banach 
space X and its topological dual X* , the sequential Painleve-Kuratowski upper/outer limit 
as x — > x is defined by 

LimsupG(x) := <x* £ X* 3 sequences Xk — > x, x* k H> x* as k — > oo 

1 " . (1.5) 

such that x* k eG{x k ) for all k € M := {1, 2, . . .}|, 

where w* signifies the weak* topology of X*. For a set SI C X the symbol x x means 
that x —> x with x € £1. If ip: X — > M := (—00,00] is an extended-real-valued function 
finite at x, the symbol x4i signifies the convergence x — > x with (p(x) — > <p(x). 



2 Tools of Variational Analysis 

In this section we briefly review some basic constructions and results of the generalized 
differentiation theory in variational analysis that are widely used in what follows. The reader 
can find all the proofs, discussions, and additional material in the books [6, 13, 14, 18, 19] 
and the references therein in both finite and infinite dimensions. Unless otherwise stated, 
all the spaces under consideration are Banach with the norm || • || and the canonical pairing 
(■, ■) between the space in question and its topological dual. 

Let us start with convex functions ip: X — > M. Given x 6 dom<^ := {x € X\ ip(x) < 00}, 
the subdifferential (collection of subgradients) of tp at x in the sense of convex analysis is 

dip(x) := {x* G X*\ {x*,x-x} < (p(x) - ip(x) for all x€X}. (2.1) 

Directly from definition (2.1) we have the following nonsmooth counterpart of the classical 
Fermat stationary rule for convex functions: 

x is a minimizer of p if and only if 6 d(p(x). (2-2) 
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The sub differential of convex analysis (2.1) satisfies a number of important calculus rules 
that are mainly based on separation theorems for convex sets. The central calculus result 
is the following Moreau-Rockafellar theorem for representing the subdifferential of sums. 

Theorem 2.1 (subdifferential sum rule for convex functions). Let <fi : X — > 1R, 

i = l,...,m, be convex lower semicontinuous functions on a Banach space X. Assume 
that there is a point x G n™ =1 dom(^j at which all (except possibly one) of the functions 
ipi, . . . ,(p m are continuous. Then we have the equality 

m m 
1=1 i=l 

Given a convex set f2 C X and a point x G J), the corresponding geometric counterpart 
of (2.1) is the normal cone to at x defined by 

N(x; fi) = jx* G X* I (x* , x — x) < for all ieQ}, (2.3) 

which is in fact the subdifferential (2.1) of the set indicator function 5(x;Q) at x that is 
equal to for x G $7 and to 00 for x ^ f2. 

Besides the aforementioned convex constructions suitable for the study of the generalized 
Fermat-Torricelli problem (1.3) in the case of convex sets fij, we need in what follows their 
extensions to nonconvex objects. 

Given an arbitrary extended-real- valued function <p: X — > M finite at x and given e > 0, 
define first the e- subdifferential of ip at x by 

d £ (f(x) := {x* G X* 

For e = the set d<p(x) := 9o ( / 9 (^) i s known as regular /viscosity /Frechet subdifferential 
of <y3 at x; it reduces to the classical gradient {Vc^(x)} when ip is Frechet differentiable 
at this point and to the subdifferential (2.1) when cp is convex. However, being naturally 
and rather simply defined, the Frechet subdifferential and its e-enlargements (2.4) do not 
possess — apart from locally convex settings and the like — a number of required calculus and 
related properties. For example, dip(x) may often be empty (e.g., for <p(x) = —\x\) and an 
analog of the sum rule from Theorem 2.1 does not hold for dcp(x) whenever e > 0; e.g., in 
the case of <pi(x) = \x\ and ^2{x) = —\x\. 

The situation dramatically changes when we employ a sequential regularization of the 
e-subdifferentials (2.4) defined by 

dip(x) := Limsup9 e v3(x) (2.5) 

■jl — rx 

via the sequential outer limit (1.5) and known as the basic /limiting /Mordukhovich subdiffer- 
ential of ip at x. We can equivalently put e = in (2.5) if (p is lower semicontinuous around 
x and if the space X is Asplund, i.e., each of its separable subspaces has a separable dual. 



lim inf -j*.*-*; > _,}. (2 . 4) 



x — x 
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The latter subclass of Banach spaces is sufficiently large including, in particular, every re- 
flexive space and every space with a separable dual. On the other hand, it does not contain 
some classical Banach spaces important for applications as, e.g., C[0, 1] and -^[0, 1]. 

A geometric counterpart of the subdifferential (2.5) is the corresponding (basic, limiting, 
Mordukhovich) normal cone to a set 0, C X at x G 0, that can be defined via the subdiffer- 
ential (2.5) of the indicator function N(x;Q) := d5(x;Q) and reduces to the normal cone 
of convex analysis (2.3) for convex sets Q. The given definition of our basic normal can be 
equivalently rewritten in the limiting form 

N(x;Q) = Limsup N £ (x;Q) (2.6) 

X— >x 

with the sets of e-normals N e (-; Q) defined for e > by 

N £ (x;Q) := {x* G X*| limsup ^ ,X _~f < e}, x G fi, (2.7) 

x^Xx 

where N £ (x;Cl) := if x ^ Q for convenience. When the set f2 is locally closed around x 
and the space X is Asplund, we can equivalently replace N £ (-;Q) in (2.7) by the prenor- 
mal/Frechet normal cone N (■•,&,) := Nq(-;Q). Furthermore, in the case of X = M n the 
normal cone (2.6) admits the representation 

N(x;n) = Limsup [cone (x -U(x;Q))], (2.8) 

X— >x 

where H(x; Q) denotes the Euclidean projection of the point x G M n onto the closed set O,, 
and where coneil signifies the collection of rays spanned on £1. Representation (2.8) was 
actually the original definition of the limiting normal cone in [12]. 

In spite of the nonconvexity of the limiting constructions (2.5) and (2.6), they enjoy 
well-developed calculus rules that are pretty comprehensive in the Asplund space setting 
and are based on variational/extremal principles; see, e.g., [13]. In particular, the following 
sum rule for the subdifferential (2.5) is used in this paper. 

Theorem 2.2 (subdifferential sum rule for nonconvex functions). Let ipi : X — > 1R, 

i = 1, . . . ,m, be lower semicontinuous functions on an Asplund space X. Suppose that all 
(except possibly one) of them are locally Lipschitzian around x G n^j^dom^j. Then we have 
the inclusion 

n n 
i=l i=l 

3 Generalized Differentiation of Minimal Time Functions 

This section is devoted to reviewing, for the reader's convenience, some recent results on 
generalized differentiation of the minimal time functions of type (1.1) developed in detail 
in our separate paper [16], which in fact was mainly motivated by the application to the 
generalized Fermat-Torricelli problem (1.3) given in what follows. Let us present and discuss 
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the major required results on generalized differentiation of the minimal time functions in 
both convex and nonconvex cases. We say that x G X is a in- set point for the minimal time 
function (1.1) if x G SI and that x is an out- of- set point for (1.1) if £ ^ SI. 

The following result, which is a consequence of Proposition 4.1, Proposition 4.2, and 
Theorem 5.2 from [16], provides precise relationships between basic subdifferential (2.5) 
of the minimal time function (1.1) and the basic normal cone (2.8) to the corresponding 
targets in the case of in-set points. 

Theorem 3.1 (basic subgradients of minimal time functions and basic normals 
to targets at in-set points). Let x G SI for the minimal time function (1.1) on a Banach 
space X , and let the support level set C* C X* be defined by 

C* := {x* G X*\ a F (-x*) < 1} (3.1) 

via the support function of the dynamics given by 

a F {x*) := sup(x*,x), x* G X*. (3.2) 

Then we have the subdifferential upper estimate 

dr£(x) c N{x;n)nc*. (3.3) 

Furthermore, the latter holds as the equality 

dT£(x) = N(x-,n)nc* (3.4) 

provided that the target set Q is convex. 

The next result gives an upper estimate of the basic subdifferential of the generally 
nonconvex minimal time function (1.1) at out-of-set points via basic subgradients of the 
corresponding Minkowski gauge 

p F {x) : = inf {t > 0| x G tF}, x G X, (3.5) 

associated with the dynamics F and via basic normals to the target set S7 at points belonging 
the the minimal time/ generalized projection of x ^ S7 to S7 defined by 

Il%(x):=(x + Tg(x)F)nn. (3.6) 

It is easy to see that the generalized projection (3.6) reduces to the standard (metric) one 
when F = ]B, i.e., when (1.1) becomes the distance function (1.2). 

To proceed, we recall that the minimal time function (1.1) is well posed at x G" $1 with 
Tq(x) < oo if for any sequence xu — > x with T^(xk) — > Tq(x) as k — > oo there is a sequence 
of projection points G n^(xfc) containing a convergent subsequence. This property is 
defined and discussed in [16]: cf. also [13, Subsection 1.3.3] for the case of distance functions. 
The following conditions are sufficient for well-posedness: 

• The target fi is a compact subset of X; 
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• The space X is finite-dimensional and Q is a closed subset of X; 

• X is reflexive, S7 C X is closed and convex, and the Minkowski gauge associated with F 
generates an equivalent Kadec norm on X, i.e., such that the weak and norm convergences 
agree on the boundary of the unit sphere of X. 

Here is the aforementioned result; cf. [16, Theorem 6.3]. 

Theorem 3.2 (basic subgradients of minimal time functions at out-of-set points 
via projections). Let x ^ $7 with Tq(x) < oo, and let the minimal time function (1.1) be 
well posed at x. Then we have the upper estimate 

dTg(x)c (J [-d PF (w-x)nN(w-n)]. (3.7) 

Finally in this section, consider the case of convexity of the minimal time function (1.1), 
which is equivalent to the convexity of its target set f2 as shown, e.g., in [16, Proposition 3.6]. 
In this case we have some specific results, which are not satisfied for general nonconvex 
minimal time functions; see [16] for more details. In particular, the convex case allows 
us to establish important connections between the basic subdifferential of (1.1) and the 
corresponding normal cone to the target enlargements 

n r := {x G X\ T$ (x) < r}, r > 0, (3.8) 

at out-of-set points x ^ $7. The following result taken from [16, Theorem 7.3] contains what 
we need for applications to the generalized Fermat-Torricelli problem in this paper. 

Theorem 3.3 (subgradients of convex minimal time functions at out-of-set points). 

Let the minimal time function (1.1) be convex, and let x ^ Q be such that Uq(x) ^ with 
r = Tq(x) < oo in (3.8). Then for any w G II^(x) we have the relationships 

dT£(x) =N(x;n r )n[-d PF (w-x)] 

C N(w; n)n[- dp F (w -x)]. { ' ' 

If in addition G F , then the inclusion in (3.9) holds as equality, and thus 

dT^ (x) = N(w; n) n [ - d PF (w - x)] . (3.10) 

4 Generalized Fermat-Torricelli Problem: Optimality Con- 
ditions in Finite and Infinite Dimensions 

This section mainly concerns qualitative aspects of the generalized Fermat-Torricelli prob- 
lem (1.3) related to deriving necessary as well as necessary and sufficient conditions for its 
solutions in convex and nonconvex cases. We also show that the obtained qualitative results 
allow us to explicitly find generalized Fermat-Torricelli points in some remarkable settings. 

Let us first establish sufficient conditions for the existence of optimal solutions to the 
generalized Fermat-Torricelli problem under consideration. 
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Proposition 4.1 (existence of optimal solutions to the generalized Fermat-Torricelli 
problem). In addition to the standing assumption of Section 1, suppose that at least one 
of the sets Q±, . . . , £l n in (1.3) is bounded and that inf xe x T(x) < oo. Then the generalized 
Fermat-Torricelli problem (1.3) admits an optimal solution in each of the following settings: 

(i) The space X is finite- dimensional. 

(ii) The space X is reflexive and all the sets Slj, i = 1, . . . ,n, are convex. 

Proof. To justify (i), suppose that the set fii is bounded. Denoting a := ini x£ x T(x), we 
immediately observe that 

{x € X\ T(x) < a + 1} C {x e X\ T^(x) < a + l}, 

and hence the level set {x £ X\ T(x) < a + 1} is bounded. By [16, Proposition 3.5] the 
minimal time function (1.1) is lower semicontinuous under the assumptions in (i). Thus we 
deduce the existence of solutions to (1.3) from the classical Weierstrass theorem. 

To proceed with the proof of (ii), recall that every convex, bounded, and closed subset of 
a reflexive space is sequentially weakly compact. Furthermore, Proposition 3.5 of [16] yields 
the lower semicontinuity of the minimal time function (1.1) with such a target set. This 
implies the weak lower semicontinuity of T(-) from (1.3) under the convexity assumptions 
made and hence ensures the existence of optimal solutions to (1.3) in case (ii) by applying 
the Weierstrass theorem in the weak topology of X to this problem. A 

It is not hard to illustrate by examples that all the assumptions made in Proposition 4.1 
are essential for the existence of optimal solutions to (1.3) . Consider for instance a particular 
case of (1.4) with X = M 2 , n = 2, Q x := {(x,y) G M 2 | y > e x }, and tt 2 := M x {0}. It is 
clear that this problem does not have an optimal solution. 

Let us further proceed with deriving optimality conditions for the generalized Fermat- 
Torricelli problem (1.3). Define the sets 

Mx)--= U [-dpF(u-x)nN(u;Qi)], xeX, (4.1) 
«eng.(x) 

provided that U^.(x) / 0. As in the proof of Proposition 4.1 (with no need of boundedness 
of the target set Q while under the standing assumptions on the dynamics F) , we can deduce 
from the generalized projection definition (3.6) that n^(x) / for all x £ U in each of the 
two following cases: 

• X is finite-dimensional and Q is closed; 

• X is reflexive, is closed and convex. 

Furthermore, it easy to observe from the construction in (4.1) that 

Ai(x) = N(x;Cl) n C* as xeQi, i = l,...,n, (4.2) 

for arbitrary closed sets fij, where the support level set C* is defined in (3.1). Useful 
relationships between the sets Ai(x) and the subdifferential dT^.(x) in the out-of-set case 
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x G Vti follow from Theorem 3.2 and Theorem 3.3 for convex and nonconvex targets fij. 
These relationships are widely used in the sequel. 

We first establish necessary optimality conditions for the general nonconvex problem 
(1.3) in infinite dimensions. For simplicity we assume that G intF, which ensures the 
Lipschitz continuity of the minimal time function (•) for all sets fij, i = l,...,n, in 

(1.3) and the possibility to apply the sum rule from Theorem 3.2. Our approach allows us 
to treat, with some elaboration, the non-Lipschitzian case when int F = by using more 
involved subdifferential formulas for the minimal time function obtained in [16] and the 
basic subdifferential sum rules for non-Lipschitzian functions given in [13, Chapter 3]. 

Theorem 4.2 (necessary optimality conditions for the generalized Fermat-Torricelli 
problem). Let X be an Asplund space, and let G int F. If x G X is a local optimal so- 
lution to the generalized Fermat-Torricelli problem (1.3) such that for each i = 1, . . . ,n the 
minimal time function Tq. (•) is well posed at x when x ^ fij, then 

n 

Oe^ii(i) (4.3) 
i=l 

with the sets Ai(x), i = 1, . . . ,n, defined in (4.1). 

Proof. If x is a local solution to (1.3), then G dT(x) by the generalized Fermat stationary 
rule; see [13, Proposition 1.114]. It is well known that the minimal time function (1.1) is 
Lipschitz continuous on X provided that G intF; see, e.g., [9, Lemma 3.2]. Employing 
thus the nonconvex subdifferential sum rule from Theorem 2.2, we have 

n 

0€ (4.4) 

i=i 

Comparing inclusion (3.3) from Theorem 3.1 with formula (4.2) gives us the upper estimate 

dT^(x)cA t (x), i = l,...,n, (4.5) 

in the in-set case x G fij. Furthermore, by Theorem 3.2 the above inclusion (4.5) holds 
also in the out-of-set case x ^ Qi under the assumed well-posedness. Substituting (4.5) into 

(4.4) , we arrive at (4.3) and complete the proof of the theorem. A 

For the particular case (1.4) of problem (1.3) in Hilbert spaces a more explicit coun- 
terpart of (4.3) holds, which provides a convenient necessary optimality condition for the 
Steiner-type extension (1.4) of the Fermat-Torricelli problem. 

Theorem 4.3 (necessary optimality conditions for the Steiner-type extension of 
the Fermat-Torricelli problem in Hilbert spaces). Let X be a Hilbert space, and let 
x G X be a local optimal solution to problem (1.4) such for each i = 1, . . . ,n the distance 
function d(-;Oj) is well posed at x when x £ £V Then condition (4.3) is necessary for 
optimality of x in (1.4), where the sets Ai(x) are explicitly expressed by: 

x — II(x; fij) 



Ai{x) 



a[X, ' ll) foralli = l,...,n. (4.6) 

N(x;ni)nlB if iefij 
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Proof. Observe that a(-) = || • || in (3.2) and that C* = IB in (3.1) by X* = X and F = B. 
By Theorem 4.2 it remains to prove that the expression for Ai(x) in (4.6) for x ^ Vti reduces 
to (4.1) in the setting under consideration. 

Fix an arbitrary vector x ^ and show that any vector u G Ai(x) from (4.1) belongs 
to the set on the right-hand side of (4.6). It is well known that in Hilbert spaces we have 

(x — oj, x — ui) < — \\x — uj\\ 2 for all u> G IT(x; f2j) and x G f2j. 

The latter implies, by definitions (2.7) and (2.6), that 

x - uj £ N(uj; Qi) C N(cu; 0^). 

Since Pf(x) = \\x\\ for the Minkowski gauge (3.5) in this case, it gives 

-dp(u;-x) = {-^-}. (4.7) 
I ||x — oj\\ J 

Using now (4.7) and the inclusion 

G N(oj; n { ) 



x — CO 



held by representation (2.8), we have for u G Ai(x) in (4.1) the relationships 

u G -dp F {u) - x) n N(uj;fli) = (-^ — ^-j. 

I \\x — uj\\ ) 

It follows that \\x — lo\\ = d(x; Qi) due to ui G II(x; Oj), and thus 



x — u) x — II(x; Oj) 

U = G 



\x — oj\\ d(x; 0;) ' 
which justifies that u belongs to the set on the right-hand side of (4.6). 

x — n(ar 

To prove the converse inclusion, take any u G _ , ' , 1 and find ui G II(x; Qj) such 

d(x; ili) 

that u = -p — p| . Then u G iV(a;; O^) by (2.8) and u G -dp F {u) - x) by (4.7). This gives 

u G —8pf{oj — x) n iV(o;; J7j) 

and shows that the set Ai(x) in (4.6) belongs to the one in (4.1) for each i G {1, . . . ,n}, 
which thus completes the proof of the theorem. A 

Observe that the well-posedness assumption of Theorem 4.3 is automatic when either 
X is finite-dimensional or the corresponding set Ui is convex. Furthermore, we also have 
II(a;; J7;) ^ in the same settings. 

Next we employ Theorem 4.3 to specify optimal solutions to (1.4) with n = 3 therein. 
Note that the condition (u,v) < —1/2 obtained in what follows means that the angle 
between these two vectors is larger than or equal to 120°, which is the crucial case in the 
classical Fermat-Torricelli problem. 
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Corollary 4.4 (necessary conditions for the generalized Fermat-Torricelli prob- 
lem with three nonconvex sets in Hilbert spaces). Let n = 3 in the framework of 
Theorem 4.3, where Qi,Q 2 , ^3 are pairwise disjoint subsets of X. The following alternative 
holds for a local optimal solution x G X with the sets Ai{x) defined by (4.6): 

(i) The point x belongs to one of the sets f2j, say Vt\, and does not belong to the two 
others. Then there are a 2 G A 2 (x) and 03 G A 3 (x) such that 

(02,03) < -1/2 and - a 2 - a 3 G N(x; Qi). (4.8) 

(ii) TTte point x does not belong to all the three sets Q\, Q 2 , and Q 3 . Then there are 
a-i G Aj(x) as i = 1, 2, 3 suc/i that 

( ai , aj ) = -1/2 /or i + j as i,j G {1,2,3}. (4.9) 

Proof. Since the sets Qi are pairwise disjoint, the settings in (i) and (ii) fully describe the 
possible location of x. Then Theorem 4.3 ensures that 

£ A^x) + A 2 (x) + A 3 (x) (4.10) 

with the sets Ai(x) defined by (4.6). Considering first case (i), we get by (4.6) and (4.10) 
vectors 02 G A 2 (x) and 03 G A${x) satisfying the relationships 

IKH = || 0.3 1| = 1 and - o 2 - a 3 G N(x;Qi) n B. (4.11) 

Due to the obvious identities 

||«2 + «3 II 2 = \\a2W 2 + 2(a 2 ,a 3 ) + ||a 3 || 2 = 2 + 2(a 2 ,a 3 ), 

the condition ||a2 + a.3 1| G IB <^=^ ||a2 + a.3 1| 2 < 1 is equivalent to (02,03) < —1/2. Thus 
the necessary optimality condition (4.10) can be equivalently rewritten in form (4.8), which 
completes the proof in case (i). 

Considering next case (ii), we get by (4.10) and (4.6) vectors a, G Ai(x) for i = 1,2,3 
satisfying the relationships 

Il a i|| = H02II = II03II = 1 an d 01+02 + 03 = 0. 

The latter implies that a\ + 02 = — 03, and hence 

(ai,o 3 ) + (02,03) = (ai + a 2 ,a 3 ) = -(a 3 ,a 3 ) = -1. 

Similarly we arrive at the conditions 

(ai,a 2 ) + (01,03) = -1 and (ai, o 2 ) + (a 2 , a 3 ) = -1. 

Subtracting pairwisely the above equalities gives us 

(01,02) = (02,03) = (ai,o 3 ) = -1/2, 

which can be written as (4.9) and thus completes the proof of the corollary. A 

The following example illustrates the application of Corollary 4.4 to a particular problem 
on the plane with two convex and one nonconvex sets. 
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Figure 1: A Nonconvex Generalized Fermat-Torricelli Problem. 

Example 4.5 (nonconvex generalized Fermat-Torricelli problem on the plane). 

In the setting of Corollary 4.4, let fii be the ball centered at c\ := (0, —2) with radius r = 1, 
let Q 2 be the ball centered at c 2 := (0, —6) with the same radius r = 1, and let £1% be a 
nonconvex set defined by 

Q 3 := {(x 1 ,x 2 ) G M 2 \ x 2 > -|xi|} (4.12) 

as depicted on Figure 1. By Proposition 4.1 there is an optimal solution to this problem. 
Applying then Corollary 4.4, all the assumptions of which are satisfied, we find two points 
lying on the boundary of J7i and denoted by u and v such that \c\ , u] is the angle bisector 
for the angle formed by the lines uc 2 and up u , where p u is the projection of u to O3, while 
[v,c\] is the angle bisector for the angle formed by the lines vc 2 and vp v . These two points 
satisfy (i) of Corollary 4.4 and they are actually the optimal solutions to the problem under 
consideration. It is not hard to find u and v numerically; we get u = (—0.8706, —2.4920) 
and v = (0.8706, —2.4920) up to five significant digits, with the optimal value of the problem 
equal to 3.7609. 

We continue now by considering the generalized Fermat-Torricelli problem (1.3) with 
convex target set Oj as i = 1, . . . , n in Banach spaces. In this case we derive necessary and 
sufficient optimality conditions for Fermat-Torricelli points. 

It follows from Theorem 3.1 and Theorem 3.3 that in the case of convex sets Qj with 
Hq. (x; fti) ^ as x ^ f2, and G F we have the equalities 

Ai(x) = -dp F (uj - x) n N(cj;Ui) for any x£X and u e U^,(x; Ui) (4.13) 

for the sets Ai(x) defined in (4.1), where the subdifferential and normal cone are explicitly 
computed by formulas (2.1) and (2.3) of convex analysis. Here is a characterization of 
Fermat-Torricelli points for convex problems. 

Theorem 4.6 (necessary and sufficient conditions for generalized Fermat-Torricelli 
points of convex problems in Banach spaces). Let all the target sets fij be convex, 
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let G inti 7 /or problem (1.3) formulated in a Banach space X, and let x £ X be such 
that Hq.(x) whenever x ^ Jlj as i = 1, . . . ,n. T/ien condition (4.3) iwtt f/ie seis -Aj(x) 
defined in (4.13) is necessary and sufficient for optimality of x in this problem. 

Proof. As mentioned above, the convexity of all the sets fl{ implies the convexity of the cost 
function T(x) in problem (1.3). By the generalized Fermat rule (2.2) for convex functions 
we get the inclusion G dT{x) as a necessary and sufficient conditions for optimality of 
x G X in (1.3). Since all the functions Tq, (•) are locally Lipschitzian under the inferiority 
assumption on the dynamics F, the convex subdifferential sum rule of Theorem 2.1 ensures 
that the latter inclusion is equivalent to 

n 

0eY,9T£(x). (4.14) 

i=i 

Applying now relationship (4.2) and equality (3.4) of Theorem 3.1 in the in-set case as well 
as Theorem 3.3 in the out-of-set case, we conclude that 8Tq.(x) = Ai(x) as i = 1, . . . , n. 
Thus inclusion (4.14) is equivalent to (4.3), and the latter is necessary and sufficient for 
optimality of x in the convex Fermat-Torricelli problem (1.3). A 

The following consequence of Theorem 4.6 provides an explicit characterization of Fermat- 
Torricelli points in the convex Steiner-type extension (1.4) of the classical problem in the 
Hilbert space setting. In this case we use formula (4.6) for constructing the sets Ai(x), 
which reduce to singletons if x £ fij and are computed explicitly by (2.3) if x G fij. 

Corollary 4.7 (characterization of optimal solutions to the convex Steiner-type 
extension of the Fermat-Torricelli problem in Hilbert spaces). Let X be a Hilbert 
space, and let all the sets in (1.4) be convex. Then condition (4.3) with A,i{x) computed 
in (4.6) is necessary and sufficient for optimality of x G X in problem (1.4). 

Proof. It follows from Theorem 4.6 due the fact that the sets A^x) from (4.1) reduce to 
those in (4.6) as proved in Theorem 4.3 and due to the projection nonemptiness II(x; fij) / 
for any x G X and i = 1, . . . , n in the setting under consideration. A 

Note that, in contrast to problem (1.4) addressed in Corollary 4.7, the characterization 
of generalized Fermat-Torricelli points for problem (1.3) obtained in Theorem 4.6 depends 
on the dynamics F and clearly determines different solutions for the same targets sets 0^ 
while different dynamics sets F. For example, consider the case of the three singletons 
fti : = {(-1,0)}, ft 2 := {(0,1)}, and fl 3 := {(1,0)} on the plane M 2 . Then Corollary 4.7 
gives us the unique optimal solution (0, l/VS) to the corresponding problem (1.4), while for 
F = [—1,1] x [—1,1] and the same sets f2j, i = 1, 2, 3, we have the unique optimal solution 
(0, 1) to the generalized Fermat-Torricelli problem (1.3) with these sets F and fij. 

Let us present a simple application of Corollary 4.7 to a version of the generalized 
Fermat-Torricelli problem (1.4) for finitely many disjoint closed intervals of the real line. 

Proposition 4.8 (Fermat-Torricelli problem for closed intervals of the real line). 

Consider problem (1.4) with the sets fij given by n disjoint closed intervals [a^fej] C M as 
i = 1, . . . , n, where a\ < b± < a<i < hi < . . . < a n < b n . The following hold: 
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(i) If n = 2k + 1, then any point of the interval [0^+1,6^+1] (the mid interval) is an 
optimal solution to the problem under consideration. 

(ii) If n = 2k, then any point of the interval [bk,a,k+i] is an optimal solution to the 
problem under consideration. 

Proof. Let <p(x) := d(x;Q) with ft = [a,b]. It is easy to compute (as, e.g., a particular 
case of Theorem 3.3) the subdifferential of 92 by: 

{0} if a < x < b, 

[-1,0] if x = a, 

[0,1] if x = b, 

{-1} if x < a, 

{1} if x > b. 

Consider first case (i) when n = 2k + 1, we get for any x G [0^+1,6^+1] the relationships 

k 

Y J dd(x;n i ) = {k}, 
i=i 

n 

dd(x;Q i ) = {-k}, 

i=k+2 

and G dd(x; flk+i)- 

The latter implies that G Yli=i dd(x; fij), which ensures by Corollary 4.7 that x is an 
optimal solution to the problem under consideration. Taking further any x ^ [0^+1,6^+1], 
we get by the above calculation that ^ ^27=1 9d(x; fij) and hence learn from the char- 
acterization of Corollary 4.7 that such a number x cannot be an optimal solution to the 
problem. The even case of n in (ii) is treated similarly. A 

Another application of Corollary 4.7 provides complete characterizations of Fermat- 
Torricelli points for the convex problem (1.4) with n = 3 in Hilbert spaces. Note that in 
this case, due the projection uniqueness, we have 

for the sets Ai(x) defined in (4.6). 

Proposition 4.9 (characterizations of generalized Fermat-Torricelli points for 
three convex sets in Hilbert spaces). Let X be a Hilbert space, and let ^1,^2, ^3 
be pairwise disjoint convex subsets of X. Then x G X is an optimal solution to problem 
(1.4) generated by these sets if and only if one of the conditions (i) and (ii) of Corollary 4.4 
is satisfied, where the vectors ai, i = 1,2,3, are defined in (4.15), and where the normal 
cone N(x;Qi) in (4.8) is computed by (2.3). 

Proof. The necessity part of the proposition follows from Corollary 4.4 by the observations 
that the convex problems under consideration is well posed at x and that 7^ fl(x; fij) is 



d<p(x) = < 
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a singleton for any i = 1,2,3. The sufficiency part of the proposition can be derived from 
Corollary 4.7 by the arguments developed in the proof of Corollary 4.4. A 

Finally in this section, we illustrate the application of Proposition 4.9 to some particular 
problems of Fermat-Torricelli type formulated on the plane. 




Figure 2: A Convex Generalized Fermat-Torricelli Problem. 



Example 4.10 (convex generalized Fermat-Torricelli problems on the plane). Let 

the sets Oi, £2 2 , and ^3 in problem (1.4) are closed balls in M 2 of radius r = 1 centered at 
the points (0,2), (—2,0), and (2,0), respectively; see Figure 2. We can easily see that the 
point (0, 1) G f^i satisfies all the conditions in Proposition 4.9(i), and hence it is an optimal 
solution (in fact a unique one) to this problem. 

More generally, consider problem (1.4) in M 2 generated by three arbitrary pairwise 
disjoint disks denoted by f2j, i = 1, 2, 3. Let ci, c 2 , and C3 be the centers of the disks. Assume 
first that either the line segment [c 2 , C3] G M 2 intersects Qi, or \c\, C3] intersects f^, or [c\, c 2 ] 
intersects O3. It is not hard to check that any point of the intersections (say of the sets f2i 
and [02, C3] for definiteness) is an optimal solution to the problem under consideration, since 
it satisfies the necessary and sufficient optimality conditions of Proposition 4.9(i). Indeed, 
if x is such a point, then a 2 and 03 from (4.15) are unit vectors with (02,03) = —1 and 
-a 2 - a 3 = G N(x; Qi). 

If the above intersection assumptions are violated, we define three points qi,q 2 , and 
qs as follows. Let u and v be the intersections of [01,02] and [01,03] with the boundary 
of the disk centered in c\. Then we can see that there is a unique point q\ on the minor 
curve generated by u and v such that the measures of angle c\q\c 2 and c\q\c^ are equal. 
The points q 2 and q% are defined similarly. Proposition 4.9 yields that whenever the angle 
C2Q1C3, or ci(?2C3, or c 2 q^c\ equals or exceeds 120° (say the angle c 2 q\c-$ does), then the point 
x := qi is an optimal solution to the problem under consideration. Indeed, in this case a 2 
and 03 from (4.15) are unit vectors with (a 2 , 03) < —1/2 and —a 2 — 03 G N(x; Qi) because 
the vector —a 2 — 03 is orthogonal to fix. 

If none of these angles equals or exceeds 120°, there is a point q not belonging to fij 
as % = 1,2,3 such that the angles c\qc 2 = c 2 qc^ = c^qc\ are of 120°, and q is an optimal 
solution to the problem. Observe that in this case the point q is also a unique optimal 
solution to the classical Fermat-Torricelli problem determined by the points c±,c 2 , and C3. 
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5 Generalized Fermat-Torricelli Problem in Convex Settings: 
Numerical Aspects 

The concluding section of the paper is devoted to some numerical aspects of solving the 
generalized Fermat-Torricelli problem (1.3) and its concretizations for the case of n convex 
target sets in finite-dimensional spaces. Based on the subgradient method in convex opti- 
mization and the subdifferential calculus results discussed in Sections 2 and 3, we develop 
a first-order algorithm of solving a general convex problem (1.3) and present some of its 
specifications and implementations. 

Theorem 5.1 (subgradient algorithm for the generalized Fermat-Torricelli prob- 
lem). Let Oj, i = 1, . . . , n, be convex subsets of a finite- dimensional Euclidean space X, let 
G int F, and let S ^ be the set of optimal solutions to problem (1.3). Picking a sequence 
{a k } as k G FV of positive numbers and a starting point x\ G X, consider the algorithm 

n 

Xk+i = x k - a k ^2v ik , A; = 1,2,..., (5.1) 
i=i 

with an arbitrary choice of vectors 

v ik G -dp(uj ik - Xk) n N(u ik ; tti) for some u ik G U-.(x k ) if x k fi* (5.2) 
and Vik = otherwise. Assume that 

oo oo 

a k = oo and I 2 := a| < 0. (5.3) 

k=l k=l 

Then the iterative sequence {xk} in (5.2) converges to an optimal solution for problem (1.3) 
and the value sequence 

V k :=mm{T( Xj )\j = l,...,k} (5.4) 
converges to the optimal value V in this problem. Furthermore, we have the estimate 

V k-V< T^k ' 

where < L < oo is a Lipschitz constant of the function T(-) from (1.3) on X. 

Proof. As mentioned above, the value function T(-) in (1.3) is convex and globally Lips- 
chitzian on X. By Theorem 3.1 and Theorem 3.3 the convex subdifferential of the minimal 
time functions (1.1) at x k is computed by 

N(x k ;Qi) fl {v G X\ cr F (-v) < l} if x k G fy, 

(5.5) 



dTi(xk) 



N{Lo ik ; Sli) n [ - dp(uj ik - x k )] if x k Q i: 



where C0i k G 11^. (x k ) is any generalized projection vector, i G {l,...,n}, and k G IN. 
Recalling now the subgradient algorithm for minimizing the convex function T(-), we have 

n 

x k+ i = x k - a k ^ Vk with v k G dT(x k ), 1,2,.... (5.6) 

i=i 
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The convex subdifferential sum rule of Theorem 2.1 provides the representation 



Vk 



= ^2v ik with v ik G dT^(x k ) 



i=i 



of the subgradient v k in (5.6). Substituting the latter into (5.6) gives us algorithm (5.1) 
with Vit satisfying (5.2). Employing now the well-known results on the subgradient method 
for convex functions in the so-called "square summable but not summable case" (see, e.g., 
[2]), we arrive at the conclusions of the theorem under the conditions in (5.3). A 

Note that, using the above arguments, we can similarly apply to the generalized Fermat- 
Torricelli problem the subgradient method for convex optimization in the other cases con- 
sidered in [2] with the corresponding replacements of the convergence conditions (5.3). 

Let us present a useful consequence of Theorem 5.1 in the setting of (1.3) when the 
Minkowski gauge (3.5) is differentiable everywhere but the origin; this holds, e.g., for the 
distance function (1.2). In the case under consideration we denote 



g F (x) := < 



Vp F (x) if 

(5.7) 

if x = 0. 



Corollary 5.2 (subgradient algorithm under smoothness assumptions). In the 

setting of Theorem 5.1, assume in addition that the Minkowski gauge p F (-) is differentiable 
at every point X \ {0}. Picking a sequence of positive numbers {a k } satisfying conditions 
(5.3) and given a starting point x\ G X , form the algorithm 

n 

x k+ i = x k + a k ^g F {uj k i-x k ), (5.8) 
i=i 

where 0Ji k G LT^. (x k ) is an arbitrary projection vector. Then all the conclusions of Theo- 
rem 5.1 hold true for algorithm (5.8). 

Proof. Fix i G {1, . . . ,n} and k G IN. When x k ^ fij we have C0i k — x k ^ for any 
ui ik G IT^.(x fe ). Hence dp F (w ik - x k ) = V ' p F {x k - io ik ) = g F {u ik - x k ) by (5.7), and the 
intersection —dp F {uo k — Xi k ) n N{uii k ; fij) reduces to the singleton {—g F {ooi k — x k )}, which 
we take for vi k in Theorem 5.1 when x k ^ In the other hand, for x k G £2j we get 

v ik := = -g F (0) = -g F {uj ik - x k ). 

Thus algorithm (5.8) in both cases agrees with (5.1) under the assumptions made. A 
We have a further specification of algorithm (5.8) for the convex problem (1.4). 

Corollary 5.3 (subgradient algorithm for convex Steiner-type extensions). Con- 
sider problem (1.4) with convex sets Oj, i = l,...,n, in a finite- dimensional Euclidean 
space X. Given a sequence {a k } of positive numbers satisfying (5.3) and a starting point 
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x\ G X , form algorithm (5.8) with <7f(0 computed by 
g F (uki ~ x k ) 



if Xk (£. Qi, 



if x k e fij. 

T/ten a// i/te conclusions of Theorem 5.1 are satisfied for this algorithm. 

x 

Proof. Follows from Corollary 5.2 with Pf(x) = \\x\\ and V pf(x) = - — - if x ^ 0. 



(5.9) 



A 



Now we consider some examples of implementing the above subgradient algorithms to 
the numerical solution of particular versions of the generalized Fermat-Torricelli problem. 

Example 5.4 (Fermat-Torricelli problem for disks). Consider the Steiner-type ex- 
tension (1.4) of the Fermat-Torricelli problem for pairwise disjoint circular disks in M 2 . Let 
Ci = (a,i, bi) and r^, i = 1, . . . , n, be the centers and the radii of the disks under consideration. 
The subgradient algorithm of Corollary 5.3 is written in this case as 



x k+ i =x k - ak^qik, 
i=i 



(5.10) 



where the quantities qik are given by 

' 

Qik 



if 



\x k 



< n, 



Xk 



if \\xk - Ci\\ > n. 



The corresponding quantities Vk are evaluated by formula (5.4) with 

n 

T ( X i)= Yl {W X j ~ C i\\- r i)- 

Writing a MATLAB program, we can compute by the above expressions the values of 
Xk and Vk for any number of disks and iterations. This allows us, in particular, to examine 
the convergence of the algorithm in various settings. The following table shows the results 
from the implementation of the above algorithm for three circles with centers (—2, 0), (0, 2), 
and (2,0) and with the same radius r = 1. The presented calculations are performed for 
the sequence ak = 1/k satisfying (5.3) and the starting point x\ = (5,7). 



MATLAB RESULT 


k 


Xk 


v k 


10 


(0.6224, 1.1995) 


2.7243 


100 


(0.0552, 0.9984) 


2.4741 


1,000 


(0.0047, 0.9995) 


2.4721 


10,000 


(0.0004, 0.9999) 


2.4721 


100,000 


(0.0000, 1.0000) 


2.4721 


1,000,000 


(0.0000, 1.0000) 


2.4721 
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Observe that the numerical results obtained in this case are consistent with the theoretical 
ones given in Proposition 4.9. 

For four disks centered at (0, 0), (2, 2), (1, 0), and (2, —2) and the same radius r = 1/4, 
the MATLAB program gives us the optimal point (0.8453, —0.0000) and the optimal value 
4.7141. For five disks centered at (-1, 0), (-1, 1), (0, 2), (1, 1), and (1, 0) with radius r = 1/2, 
we get the optimal solution (0.0000, 0.8505) and the optimal value equal to 3.2973. 

Next we apply the subgradient algorithm (5.8) to the Steiner-type extension (1.4) of the 
Fermat-Torricelli problem for the case of squares which is significantly different from 
the case of disks in Example 5.4. 

Example 5.5 (Fermat-Torricelli problem for squares). Consider problem (1.4) gen- 
erated by pairwise disjoint squares f2j, i = 1, . . . , n, of right position in IR 2 , i.e., such that 
the sides of each square are parallel to the x-axis or the y-axis; see Figure 3. The center 
of a square is the intersection of its two diagonals, and its radius equals one half of the 
side. Let q = (a^, 6«) and ri, i = 1, . . . , n, be the centers and the radii of the squares under 
considerations. Then the vertices of the ith square are denoted by vu = (a« + r^, hi + rj), 
V2i = (ai ~ ri, hi + ri), v 3i = (ai - ri, hi - r^), and v 4i = (a« + r i: hi - r,j). 

Given a starting point x\ and a sequence {o^} satisfying the conditions in (5.3), the 
subgradient algorithm of Corollary 5.3 can be written in form (5.10), where xu = {x\k,x 2 k) 
and where the quantities qik are computed as follows: 

if \xik - di\ < ri and \x 2 k - h\ < ri, 
if xik - di> r,i and x 2 k - h > r, h 

if xik - ^ < and x 2 k - h > ri, 

if xik - ai < -ri and x 2 k — h < -ri, 

if x lk - ai> r,i and x 2k - h < -ri, 

if \xik - (h\ < ri and x 2 k — h > n, 
if \xik - ai\ < ri and x 2k - h < -ri, 
if xik - aj > ^ and \x 2k - h\ < ri, 

if ~ ai < -ri and \x 2k - bi\ < r» 

for alH G {1, . . . , n} and k G IV, with the corresponding value sequence defined by (5.4). 

Considering the implementation of the above algorithm for three squares with centers 
(—2, 0), (0, 2), and (2, 0) and with the same radius r = 1/2, we arrive at the optimal solution 



r o 



Xk 


- vu 


\\Xk 


- V li\\ 


Xk 


- v 2i 


\\Xk 


~ V 2i\\ 


Xk 


~ V3i 


\\Xk 


~ v 3i\\ 


Xk 


- v 3i 



Xk - V4_i 



(0,1) 

(o,-i) 

(1,0) 
, (-1,0) 
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(0, 1.3660) and the optimal value 3.5981. At the same time applying the theoretical results 

\/3 ~t~ 1 

of Corollary 5.3 to this case gives us the exact optimal solution (0, — - — ) with the optimal 
2 ~\~ 3\^3 

value , which are consistent with the numerical calculation. 




Figure 3: A Generalized Fermat-Torricelli Problem for Five Squares. 

For five squares with centers at (— 1, 0), (— 1, 1), (0, 2), (1, 1), and (1,0) and the same 
radius r = 1/4 we get by the subgradient algorithm the optimal solution (0.0000,0.7242) 
with the optimal value equal to 4.3014. However, it does not seem to be an easy exercise 
to solve the above problem theoretically by using Corollary 5.3. 

Let us finally illustrate applications of the subgradient algorithm of Theorem 5.1 to 
solving the generalized Fermat-Torricelli problem (1.3) formulated via the minimal time 
function (1.1) with non-ball dynamics. For definiteness we consider the dynamics F given 
by the square [—1,1] x [—1,1] on the plane. In this case the corresponding Minkowski gauge 
(3.5) is given by the formula 

Pf{%u x 2) = max{|xi|, \x 2 \}. (5-11) 

Note that the function pf(') fads to be differentiable at every nonzero point of M 2 , so we 
have to relay on the subgradient algorithm of Theorem 5.1 but not of its corollaries above. 
Observe also that to implement algorithm (5.1) we need to know just one element from 
the set on the right-hand side of (5.2) for each % G {1, . . . , n} and k 6 JN. By Theorem 3.3 
the latter set agrees with the subdifferential of the minimal time function Tq. (x^). 

In the following proposition we compute a subgradient of the minimal time function (1.1) 
generated by the Minkowski gauge (5.11) and a square target in M 2 , which is used then to 
construct a subgradient algorithm to solve the corresponding Fermat-Torricelli problem. 

Proposition 5.6 (subgradients of minimal time functions with square dynamics 
and targets). Let F = [—1,1] x [—1,1], and let Q be a square of right position in M 2 
centered at c = (a,b) with radius r > 0. Then a subgradient v(x\,X2) G dT^{x\,X2)(not 
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necessarily uniquely defined) of the minimal time function Tq{x\,X2) at (xi,X2) is computed 
by 

(1,0) if \x2 — b\ < x\ — a, x\ > a + r, 



v(xi,x 2 ) = < 



(—1, 0) if \x2 — b\ < a — xi, x\ < a — r. 
(0, 1) if \x\ — a\ < x 2 — b, X2 > b + r, 
(0, —1) if \x\ — a\ < b — X2, X2 < b — r, 



(5.12) 



if (x±, X2) £ fi. 

Proof. It is ease to see that the minimal time function (1.1) admits the representation 

Tl(x) = info F (x;fi). 

This implies by (5.11) and the above structure of fl that 

X\ — (a + r) if \x 2 — b\ < x\ — a, x\ > a + r, 

(a — r) — x\ if |x2 — 6| < a — xi, xi < a — r, 

^2 — (b + r) if \x± — a\ < X2 — b, X2> b + r, 

{b — r) — X2 if \x\ — a\ < b — x 2 , X2 < b — r, 

if (xi,x 2 ) G ^ 

for all (x±,x 2 ) G -K 2 . Applying to (5.11) the well-known subdifferential formula for maxi- 
mum functions in convex analysis allows us to compute 



Tg( Xl ,X2) = < 



dp F (xi,x 2 ) = < 



{v u v 2 ) G M 2 


M + 


M 


<!} 


if (x u x 2 ) = (0,0), 


(0,1)} 








if 


\xi\ < x 2 , 


(o,-i)} 








if 


x 2 < —\xi\, 


(1,0)} 








if 


xi > \x 2 \, 


(-1,0)} 








if 


xi < ~\x 2 \, 


{v u v 2 ) G 


M + 


M 


= 1, V! > 0, w 2 > 0} 


if 


x\ — x 2 > 0, 


(«1,W2) G 


M + 


M 


= 1, Vl > 0, v 2 < 0} 


if 


xi = -x 2 > 0, 


(«1,W2) G J? 2 


M + 


M 


= 1, Vl < 0, w 2 < 0} 


if 


x\ — x 2 < 0, 


(«1,W2) G J? 2 


M + 


M 


= 1, Vl < 0, v 2 > 0} 


if 


xi = -x 2 < 0. 
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In this way we can check by Theorem 3.3 that the vector v(xi,x 2 ) is a subgradient of 
T(-, •) at (x, x 2 ), which completes the proof of the proposition. A 

Now we are able to implement the subgradient algorithm of Theorem 5.1 to the problem 
under consideration. 

Example 5.7 (implementation of the subgradient algorithm). Consider the gen- 
eralized Fermat-Torricelli problem (1.3) with the dynamics F = [—1,1] x [—1,1] and the 
square targets fij of right position centered at (<2j, 6j) with radii T{ as i = 1, . . . , n. Given 
a sequence of positive numbers {a k } satisfying (5.3) and a starting point xi, construct the 
subgradient algorithm (5.1) for the iterations Xk = {xik,x 2k ) in Theorem 5.1, where the 
vectors are computed by Proposition 5.6 as 

(1, 0) if \x 2 k ~ bi\ < xik - CH and x lk > a; + n, 
(-1,0) if \x 2 k ~ bi\ < a,i - xik and Xi k < a* - n, 
< (0, 1) if \xik - cij| < x 2 k - h and x 2 k > h + r i: 
(0,-1) if \x ik - a, I < bi - x 2 k and x 2k <bi-ri, 
(0, 0) otherwise. 

Implementing this algorithm for the case of three squares centered at (—2, 0), (0, 2), and 
(2,0) with radius r = 1/2, = 1/k and the initial point (1, 1), we arrive at an optimal 
solution (0.0000, 1.5000) and the optimal value equal to 3.0000. For five squares centered 
at (—1,0), (—1, 1), (0,2), (1, 1), and (1,0) with radius r = 1/4, we have the optimal solution 
(0.0000, 1.0000) and the optimal value equal to 3.7500. 

Acknowledgements. The authors are grateful to Jon Borwein, Michael Overton, and 
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